Angular Position of Nodes in the Superconducting Gap of Quasi-2D Heavy-Fermion 

Superconductor CeCoIns 
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The thermal conductivity of the heavy-fermion superconductor CeCoIns has been studied in a 
magnetic field rotating within the 2D planes. A clear fourfold symmetry of the thermal conduc- 
tivity which is characteristic of a superconducting gap with nodes along the (±7r, ±7r)-directions is 
resolved. The thermal conductivity measurement also reveals a first order transition at H C 2, indi- 
cating a Pauli limited superconducting state. These results indicate that the symmetry most likely 
belongs to d x 2_ y 2, implying that the anisotropic antiferromagnetic fluctuation is relevant to the 
superconductivity. 
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The superconductivity with unconventional pairing 
symmetry has been a central subject of the physics of 
superconductors. In the last two decades unconven- 
tional superconductivity has been found in several heavy 
fermion materials, such as CeCu 2 Si2, UPt 3 , UPd2Al 3 , 
UBei3 Q and the recently discovered UGe2 (|]. There, 
the relationship between magnetism and superconduc- 
tivity is the most important theme of research, because 
the strong electron-electron correlation effect which orig- 
inates from the magnetic interaction between 4/ or 5/ 
moment and itinerant electrons allows a non-phonon 
mediated pairing with unconventional symmetry. Very 
recently it has been reported that the family CeTIns 
(T=Rh, Ir, and Co) are heavy fermion superconductors 
HQ. Especially, both Celrlns and CeCoIns are ambient 
pressure superconductors, with transition temperature of 
0.4 K and 2.3 K, respectively. Subsequent observations 
of power law temperature dependence of the specific heat 
[3, 5 1, thermal conductivity Q, and NMR relaxation rate 
[7,3 have identified CeTIns as unconventional supercon- 
ductors with line nodes. The unique feature of these 
materials is that they bear some analogy with high T c 
cuprates. For example, the superconductivity appears in 
the neighbor of the antiferromagnetic state Q . Moreover, 
the crystal structure is tetragonal which can be viewed 
as layers of Celn 3 separated by layers of TIn 2 and elec- 
tronic structure is quasi-2D |^-^] . Therefore they present 
an uncommon opportunity to study the unconventional 
superconductivity in the heavy fermion materials. 

Unconventional superconductivity is characterized by 
the superconducting gap structure which has nodes along 
certain directions. Since the superconducting gap func- 
tion is intimately related to the pairing interaction, its 
identification is crucial for understanding the pairing 
mechanism. However, the detailed structure of the gap 
function, especially the direction of the nodes, is an un- 
resolved issue in most of the unconventional supercon- 
ductors. In fact, to the best of our knowledge, it is only 
in high-Tc cuprates and the B-phase of UPt3 in which 



the nodal direction has been successfully specified. The 
main reason for this is that the standard techniques to 
probe the unconventional superconductivity, such as pen- 
etration depth, specific heat, ultrasonic attenuation, and 
NMR relaxation rate do not provide direct information 
on the node directions. The most definitive test is a 
phase sensitive experiment which has been done for high- 
Tc cuprates ||. However, this technique appears to be 
available only for high-T c cuprates up to now. The situ- 
ation therefore strongly confronts us with the need for a 
powerful directional probe. 

Recently it has been demonstrated both experimen- 
tally and theoretically that the thermal conductivity k, 
which responds to the unpaired quasiparticles (QPs) be- 
low T c , is a powerful tool for probing the anisotropic gap 
structure JToj — [15^1 . An important advantage of the ther- 
mal conductivity is that it is indeed a directional probe, 
sensitive to the relative orientation among the thermal 
flow, the magnetic field, and nodal directions of the or- 
der parameter. In fact, a clear 4-fold modulation of k 
with an in-plane magnetic field which reflects the angu- 
lar position of nodes of d x 2_ v 2 symmetry has been ob- 
served in YBa2Cu3C>7_,5 [ pi3|Jll|Jl3]| , demonstrating that 
the thermal conductivity can be a relevant probe of the 
superconducting gap structure. In this Letter, we have 
measured the thermal conductivity of CeCoIn5 in mag- 
netic field rotating within the 2D CeRi3 planes. Thanks 
to the two dimensionality, a 4-fold symmetry of k which 
clearly demonstrates the presence of gap nodes in the 
(±7r, ±7r)-directions is observed. We have also found the 
first order phase transition (FOPT) at H C 2 for the first 
time. On the basis of these findings, we discuss the na- 
ture of the superconducting gap function of CeCoIns . 

Single crystal CeCoIns (T c =2.3 K) was grown by the 
self-flux method . The residual resistivity ratio (RRR) 
was approximately 18. The thermal conductivity was 
measured by the steady-stated method with one heater 
and two RuC>2 thermometers. The sample was cut into a 
rectangular shape (3.80 x 0.38 x 0.12mm 3 ) and the heat 
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current q was applied along the [100] direction. In the 
present measurements, it is very important to rotate H 
within the 2D Celn3 planes with high accuracy because 
a slight field-misalignment produces a large effect on k 
due to the two dimensionality. For this purpose, we used 
a system with two superconducting magnets generating 
H in two mutually orthogonal directions and a 3 He cryo- 
stat equipped on a mechanical rotating stage with a min- 
imum step of 1/500 degree at the top of the Dewar [ fL2"| . 
Computer-controlling two magnets and rotating stage, 
we were able to rotate H continuously within the Celn3 
planes with a misalignment of less than 0.02 degree from 
the plane, which we confirmed by the simultaneous mea- 
surement of the resistivity p. 

We first discuss the T- and H- dependences of k. The 
inset of fig. 1(a) shows the T-dependence of k and p. 
Upon entering the superconducting state, k exhibits a 
sharp kink and rises to the maximum value at T ~ 
1.7 K. The upturn of k is reminiscent of high- T c cuprates 

f. Similar T-dependence of k has been reported in Ref. 
Compared to Ref. ||, our k at the onset is slightly 
larger, but the enhancement below T c is smaller. The 
Wiedemann- Franz ratio L = ^p ~ 1.02£o at T c is very 
close to the Lorenz number L = 2.44 x 10~ 8 flW/K, in- 
dicating that the electronic contribution well dominates 
over the phonon contribution. Therefore, the enhance- 
ment of K below T c is due to the suppression of the QP 
scattering rate, similar to the high-T c cuprates. 

Figures 1(a) and (b) depict //-dependence of k for H\\ 
ab (H c2 ~ 11 T) and HL ab (H c2 ~ 5 T) below T c , 
respectively. At all temperatures, k decreases with H and 
the //-dependence becomes more gradual with decreasing 



T in both configurations. Interestingly, for HI. ab, n 
jumps to the normal state value at H c2 below 1.0 K (see 
also the inset of Fig. 1(b)). The magnitude of the jump in 
k increases with decreasing T. Since the jump in k most 
likely comes from an entropy jump, this result provides 
a strong evidence of the occurrence of a first-order phase 
transition (FOPT). As far as we know, this is the first 
material which shows a FOPT at H c2 , though a FOPT is 
predicted to occur in the Pauli paramagnetically limited 
superconducting state |l7[ ]. We will discuss this subject 
later. 

The understanding of the heat transport in the mixed 
state of superconductors with nodes has largely pro- 
gressed during the past few years H]. There, the dom- 
inant effect in a magnetic field is the Doppler shift of 
the delocalized QP energy spectrum, which occurs due 
to the presence of a supcrfluid flow around each vortex, 
and generates a nonzero QP density of states (DOS) at 
the Fermi surface fyj. While the Doppler shift increases 
K with H through the enhancement of the DOS, it can 
also lead to a decrease of K through the suppression of 
impurity scattering time and Andreev scattering time off 
the vortices. At high temperatures, the latter effect is 
predominant, while at low temperatures the former effect 
can exceed the latter effect, as demonstrated in high-T c 
cuprates pfj| . The data in Figs. 1 (a) and (b), in which 
the //-dependence of k becomes more gradual with de- 
creasing T, are consistent with the Doppler shift. Thus 
at least above 0.35 K, the Andreev scattering of the QPs 
is the main origin for the //-dependence of k. 

We now move on to the angular variation of k upon the 
rotation of H within the Celn3 planes. Figure 2 displays 
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FIG. 1. Thermal conductivity as a function of H for (a) 
H || [010] and (b) H\\ [001] below T c . The thermal current q 
is applied along [100]-direction. Inset of (a) : k and p in zero 
field. Inset of (b) : //-dependence of k near H C 2 at 0.45 K (o) 
and 0.34 K (•) 



0.65 



S 



0.60 



0.55 



1 

r=0.45K 
4// [100] 


1 1 


= 0.5T 




q-ne °\ 




oSbC _ 
SL IT 






. 2T 


V 










3 T 


1 


j H//[010] 





45 



90 
6 ( deg. ) 



135 



180 



FIG. 2. Angular variation (9 =(q,H)) of k(H, 0)/k„ at sev- 
eral H for CeCoIn 5 . The solid lines represent the result of the 
fitting by the function k(H, 6) = Co + C29 cos 29 + C^g cos 4$, 
where Co, C^e and Cab are constants. The solid circles rep- 
resent k(H, 9) at H=l T which are obtained under the field 
cooling condition at every angle. For details, see the text. 
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k(H,9) as a function of 6 =(q, H) at T=0.45 K, which 
is measured in rotating 9 after field cooling at 9 = 0°. 
The consecutive measurement inverting the rotating di- 
rection did not produce any hysteresis in k(H, 9). More- 
over, the solid circles in Fig. 2 shows k(H, 9) at if =1 T 
which are obtained under the field cooling condition at 
every angle. k(H, 9) obtained by different procedures 
of field cooling well coincide with each other. Thus the 
field trapping effect related to the vortex pinning is neg- 
ligibly small. In all data, as shown by the solid lines 
in Fig. 2, k(H, 9) can be decomposed into three terms 
with different symmetries; k{9) = kq + kiq + K^g where 
kq is a ^-independent term, and K20 = C28 cos 29 and 
Kig = C40 cos 49 are terms with 2- and 4- fold symmetry 
with respect to the in-plane rotation, respectively. The 
term K29, which has a minimum at HLq, appears as a 
result of the difference of the effective DOS for QPs trav- 
eling parallel to the vortex and for those moving in the 
perpendicular direction. 

Figures 3 (a)-(d) display K^g normalized by the normal 
state value n n after the subtraction of the kq- and K2g- 
terms from the total k. It is apparent that K^g exhibits 
a maximum at ff||[110] and [1,-1,0] at all temperatures. 
Figure 4 and the inset show the T- and H- dependences 
of \C4g\/n n . Below T c the amplitude of K^g increases 
gradually and shows a steep increase below 1 K with 
decreasing T. At low temperatures, \C^g\jK n becomes 
larger than 2%. It should be noted that this amplitude 
is more than 20 times larger than those of the 2D super- 
conductor Sr 2 Ru04 with isotropic gap in the planes [jL2j. 
Then the most important subject is "Is the observed 4- 
fold symmetry a consequence of the nodes?" . We here 
address the origin for the 4-fold symmetry. There are 
several possible origins for this. The first is the in-plane 
anisotropy of H c2 . According to Ref. ||, H c2 has very 
small but finite in-plane anisotropy; H C 2 || [100] is ap- 
proximately 2.7% larger than H C 2 || [110]. However, this 
anisotropy is too small to explain the large amplitude of 
\C4,g\/K n > 2% at H <C H C 2- Further, and more impor- 
tantly, if this 4-fold symmetry had come from the fact 
that H c2 || [100] is larger than H c2 || [110], the overall 
sign of this term should be opposite to the one actually 
observed in K^g. The second possibility is the tetragonal 
band structure inherent to the CeCoIns crystal. If the 
in-plane anisotropy of the Fermi surface is large, then the 
large anisotropy of K^g should be observed even above T c . 
However, as shown in Fig. 4, the observed 4-fold symme- 
try above T c is extremely small; \Cig\/K n <0.2 %. Thus 
the anisotropics arising from H C 2 and the band structure 
are incompatible with the data. Moreover, the amplitude 
of the 4-fold symmetry well below T c becomes more than 
10 times larger than that above T c . These considerations 
lead us to conclude that the 4-fold symmetry with large 
amplitude well below T c originates from the QP structure. 

We now address the sign of the 4-fold symmetry. In the 
presence of nodes perpendicular to the layers, the term 
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FIG. 3. (a)-(d) The 4- fold symmetry K40 / K n at several tem- 
peratures. 
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FIG. 4. Amplitude of the 4-fold symmetry \C4e\/n n as a 
function of T/T c . Inset: Same data as a function of H/H C 2- 

K^g appears as a result of two effect. The first is the 
DOS oscillation associated with the rotating H within 
the afe-plane This effect arises because the DOS 

depends sensitively on the angle between H and the di- 
rection of the nodes of the order parameter, because the 
QPs contribute to the DOS when their Doppler-shifted 
energies exceed the local energy gap. The second effect 
is the the quasiparticle lifetime from the Andreev scat- 
tering off the vortex lattice, which has the same symme- 
try as the gap function [lO|Jll]JT^| . As discussed before, 
the second effect is predominant in our temperature and 
field range. In this case, k attains the maximum value 
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when H is directed to the nodal directions and becomes 
minimum when H is directed along the antinodal di- 
rections jl(|[ll 14 1 . Thus the sign of the present 4- fold 
symmetry indicates the superconducting gap with nodes 
located along the (±tt,±tt)- directions, similar to high-T c 
cuprates. 

A quantitative comparison of the amplitude of the 
4-fold symmetry with the theory reinforces this con- 
clusion. According to Ref. ]l4j, the 4- fold symme- 
try arising from the Andreev scattering off the vor- 
tices in d-wave superconductors is roughly estimated as 

nne/nn = ~ A ^r> A cos 40. Here A is the supercon- 
ducting gap, r is the quasiparticle relaxation rate and 



'2evfv'j:H C 2/h with Vf and v'^ the in- and out- 

of-plane Fermi velocity, respectively. According to the 
numerical result, A(T) is nearly zero at T/T c >0.4 and 
shows rapid increase with decreasing T at lower tempera- 
tures. Similar tendency in the T-dependence of \Cas\/ n n 
is observed, as shown in Fig. 4. Using T ~ 1.3 x 10 11 s _1 , 
v f ~ 1 X 10 4 m/s, »/-5x 10 3 m/s @, H c2 ~ 11 T, 
2A/fc s T c =3.54, and v4(T)=0.033 at T=0.35 K from Ref. 
|L4[ gives \C4e\/nn ~8%, which is in the same order to 
the data. Thus Andreev scattering yields \Ci$\/K n which 
is consistent with the data. It is interesting to compare 
our results on CeCoIns with the corresponding results on 
YBa2Cu307_a, in which the 4-fold symmetry has been 
reported in the regime where the Andreev scattering pre- 
dominates. The observed amplitude of 4-fold symmetry 
in YBa 2 Cu 3 7 _5 is small; ~ 0.4% of total k at 6.8 K. 
However, this amplitude occupies a few % in the elec- 
tron thermal conductivity, because the phonon contribu- 
tion is about 80-90% of the total k. In YBa2Cu307_,5 
with d x 2_ y 2 symmetry, k^q has maxima at iJ||[110] and 
[1,-1,0]. Thus K4o of CeCoIn 5 is quantitatively in accord 
with ka$ of YBa2Cu307_a. 

We finally discuss the symmetry of CeCoIns inferred 
from the present results. The fact that H C 2 is determined 
by the Pauli paramagnetic limit is a direct evidence of a 
spin singlet pairing, which is consistent with the recent 
Knight-shift measurements pi] ]. Together with the fact 
that the superconducting gap has nodes at odd multi- 
ples of 45° in fc-space, we are naturally lead to conclude 
that CeCoIn$ most likely belongs to the d x 2_ y 2 symmetry 
p2fl . The d x 2_ y 2 symmetry strongly suggests that the 
anisotropic antiferromagnetic fluctuation plays an im- 
portant role for the occurrence of the superconductivity. 
This observation is in conformity with recent NMR and 
neutron scattering experiments which report anisotropic 
spin fluctuation pJZql. 

In summary, we have measured the thermal conduc- 
tivity of the quasi-2D heavy fermion superconductor 
CeCoIn.5 as a function of the relative orientation of the 
crystal axis and the magnetic field rotating within the 2D 
planes. A clear 4-fold symmetry of the thermal conduc- 
tivity which is characteristic of a superconducting gap 



with nodes at odd multiples of 45° is revealed. Rather 
surprisingly, we also observed a first order phase tran- 
sition at H C 2 at low temperatures, indicating the Pauli 
paramagnetically limited superconducting state. These 
results show that the symmetry of CeCoIn5 most likely 
belongs to d x 2_ y 2, implying that the anisotropic antifer- 
romagnetic fluctuation plays an important role for the 
superconductivity. This material is the second example 
followed by high-T c cuprates, in which the nodal struc- 
ture in the plane is successfully specified. 

Note Added. After completion of this work, we became 
aware of the work by T. Sakakibara et al. who observed 
the FOPT at H C 2 by the magnetization measurements 
II- 
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